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Extracting (11,7) direct capture cross sections from Coulomb dissociation: 

application to 14 C(n,7) 15 C 
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A methodology for extracting neutron direct capture rates from Coulomb Dissociation data is 
developed and applied to the Coulomb dissociation of 15 C on 208 Pb at 68 MeV/nucleon. Full Con- 
tinuum Discretized Coupled Channel calculations are performed and an asymptotic normalization 
coefficient is determined from a fit to the breakup data. Direct neutron capture calculations using 
the extracted asymptotic normalization coefficient provide (n, 7) cross sections consistent with di- 
rect measurements. Our results show that the Coulomb Dissociation data can be reliably used for 
extracting the cross section for 14 C(n,7) 15 C if the appropriate reaction theory is used. The result- 
ing error bars are of comparable magnitude to those from the direct measurement. This procedure 
can be used more generally to extract capture cross sections from breakup reactions whenever the 
desired capture process is fully peripheral. 

PACS numbers: 21.10.Jx, 25.60.Gc, 25.60.Tv 
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Introduction. Neutron direct capture reactions can 
play an important role in a variety of astrophysical sites, 
including big bang, stellar environments, and supernovae. 
Often, at the relevant energies, these capture rates are 
extremely small, making the direct measurement very 
difficult. Presently, direct measurements of neutron cap- 
tures on short lived unstable nuclei are not possible. The 
Coulomb dissociation method [l[ provides an alternative 
option with different systematic errors and theoretical 
challenges. The aim of this work is to provide a method- 
ical and reliable procedure to extract the desired (n, 7) 
cross section from intermediate energy Coulomb Dissoci- 
ation data. 

Here we consider the particular example of 
14 C(n,7) 15 C. This reaction is of interest to astro- 
physics for a variety of reasons: (i) it is the slowest in 
the neutron induced CNO cycle that takes places in AGB 
stars |2[; (ii) it has impact on the abundances of the 
heavy elements produced by non-homogeneous big bang 
models [3]; (iii) it modifies the abundances resulting 
from the r-process in massive Type II supernovae Q. In 
addition, it is one of the few that has been repeatedly 
studied through Coulomb Dissociation and for which 
direct measurements exist. 

Shortly after its astrophysical relevance was identified, 
a first measurement of 14 C (11,7) 15 C was performed at 
Karlsruhe @. As the container surrounding the target 
had been strongly activated by a previous experiment, 
the measurement had to be repeated many years later 
@ in order to obtain a reliable result. In the meantime, 
a number of Coulomb dissociation experiments have been 
performed using a 15 C beam on a 208 Pb target, at sev- 
eral beam energies: the lowest data is available at 35 
MeV/nucleon |8j, an intermediate energy measurement 
was performed at 68 MeV/nucleon |9J and a high energy 
experiment at 600 MeV/nucleon [XOj ] . Meanwhile, an- 
other indirect method to obtain the (n.,7) rate based on 



mirror symmetry pd| showed discrepancies between all 
three methods: the direct method, the indirect Coulomb 
dissociation method and the mirror symmetry method. 
Then recently, the new direct measurements by Reifarth 
were revised and published in a conference proceeding 
0, which agreed with the Coulomb dissociation mea- 
surements from RIKEN [9] and the mirror symmetry re- 
sults [U . The discrepancy with the MSU data [|[ still 
remained. 

Due to the weak binding of the 15 C ground state, and 
because there are no low lying resonances, the (n, 7) 
cross section is mainly determined through the El di- 
rect transition from an initial p-wave scattering state to 
the ground state, with a small branch to the only excited 
bound state [ll|. It is well known that this dominant El 
transition is completely peripheral and thus essentially 
fixed by the asymptotic normalization coefficient (ANC) 
of( 14 C| 15 C) [lj. The 1/2+ ground state is a halo n- 14 C 
s-wave, and the 5/2 + excited state is mostly d-wave. In 
both cases the relevant physics can be described by a 
single particle model. As the continuum is structureless 
below 1 MeV, the capture cross section in this low energy 
region is not very dependent on the details of the scatter- 
ing potential chosen — such as the radius and diffuseness 
- but is strongly dependent on the ANC of the bound 
state wavefunction. 

While the calculation of the neutron capture rate is 
straightforward, the Coulomb dissociation requires a reli- 
able reaction model. All three Coulomb dissociation mea- 
surements were analyzed with the virtual photon method 
[l3| . which is based on a first-order Coulomb-only semi- 
classical theory [l4j . In Q, the nuclear component was 
determined from the scaling of the breakup cross section 
on lighter targets, and subsequently subtracted from the 
total measured cross section. In [9j, the nuclear was ex- 
cluded through the selection of impact parameters larger 
than b > 30 bm. Once a Coulomb only cross section 
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is obtained, the virtual photon method provides a sim- 
ple proportionality relation to extract the desired (n, 7) 
cross section. 

The continuum discretized coupled channel method 
(CDCC) [l5| provides a non-perturbative framework in 
which to describe the breakup process, treating Coulomb 
and nuclear effects on the same footing. Multipolc ex- 
citations are fully taken into account as well as final 
state interaction effects (e.g. [Hj]). Systematic studies 
of Coulomb dissociation for loosely bound systems on a 
variety of targets, spanning a range of beam energies, 
have shown that nuclear scaling is not always reliable 
and nuclear-Coulomb interference can be very large [TtJ • 
In that work it is suggested that, rather than massaging 
the data to obtain a Coulomb only cross section, CDCC 
be used as a standard tool to analyze the full measured 
dissociation data. In this frame of mind, one would start 
with a single particle structure model for the projectile 
15 C= n+ 14 C and would adjust the parameters to ob- 
tain a good description of the breakup data. The reac- 
tion model then already takes into account nuclear and 
Coulomb effects in a coherent manner. That same poten- 
tial model for the projectile which fits the breakup data 
would then provide the corresponding neutron capture 
cross section. 

A potential model for the projectile describing the 
breakup cross sections is not unique. Indeed, as shown in 
[HI], for loosely bound systems, even if the geometries dif- 
fer significantly, it is essentially the ANC that determines 
the normalization of the breakup cross section. This is 
also true for the neutron capture cross section whenever 
it is completely peripheral, as in the case of 14 C(n,7) 15 C. 
The use of ANCs in the analysis of peripheral reactions 
was first introduced for transfer processes [Hj]. In this 
work we propose to use the dissociation data to extract 
the ANC of 15 C(g.s.) and then obtain from it the neutron 
capture cross section. We concentrate on the breakup 
data at 68 MeV/nucleon first, to show the reliability of 
our proposed procedure, and then discuss difficulties in 
the analysis of the data at 35 MeV/nucleon. We do not 
include the 600 MeV/nucleon experiment in our studies 
since, at these energies, converged CDCC calculations are 
very computational intensive and, given that relativistic 
effects are too strong to be considered approximately, re- 
sults would not be reliable. A study of the theoretical 
sources of uncertainty in the procedure will also be pre- 
sented. 

Methodology. In this work, we propose a procedure for 
extracting direct neutron capture rates for loosely bound 
systems from Coulomb dissociation, independent of the 
specifics for the experiments. It relies on the fact that 
both, the breakup cross section and the capture cross 
section are proportional to the square of the ANC. The 
method is applied to 14 C(n,7) 15 C that has been measured 
directly but also through Coulomb dissociation. 

Our starting point is a set of 14 C+n potentials which 
span a range of ANCs. In Table|T]we compile some 14 C+n 
potentials from the literature and add a few of our own, 
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TABLE I: 1 C+n potential parameters and calculated ANCs. 
The spin orbit depth was fixed at 5 MeV with the same radius 
and diffuseness as the central part. 



all with changing geometries. The V s potential strength 
is obtained from fitting the binding energy of the s-wave 
ground state of 15 C. The same V s reproduces well the 15 C 
d-wave resonance. The V p potential strength is obtained 
from the binding energy of the P3/2 neutrons in 14 C. V s 
(V p ) is used to calculate the scattering states with even 
(odd) parity. For each of these single particle potentials, 
CDCC calculations for the breakup of 15 C on 208 Pb are 
performed. Both bound states are included and fully cou- 
pled in the calculation. The CDCC model space needed 
for convergence includes partial waves for the internal 
motion of the projectile up to I < 3 and maximum 14 C+n 
relative energy E xe i < 8 MeV, with a fine discretization 
from — 2 MeV to obtain higher accuracy in this region 
(the number of bins below 2 MeV in each partial wave 
is N Sl/2 = 10, N Pl/2 = 10, N Pa/2 = 20, N d = 8,N f = 5). 
For the projectile-target relative motion, the maximum 
partial wave included is £ max = 9000. The coupling po- 
tentials are expanded in multipoles and we include all 
multipolarities up to Q — 3. As to the radial trunca- 
tions, energy bins are integrated out to 50 fm and the 
CDCC equations are solved up to i? max = 1000 fm. 

For each of the breakup energy distribution theoretical 
curves, the % 2 for the data of 0] was calculated. The 
function X 2 (C 1/2 ) is quadratic and from its minimum 
we determine the ANC Co, with an error from Xmin + 1 
[20| . Based on this range of ANCs, neutron capture cross 
sections are calculated with an error bar originating from 
the fit to the breakup data. The code FRESCO is used 
for both the breakup and capture calculations. 

In summary, the proposed procedure involves a few 
steps: (i) define a set of neutron single particle poten- 
tials that generate a range of ANCs; (ii) calculate the 
corresponding breakup energy distributions in CDCC; 
(iii) for each single particle potential, calculate the x 2 
to the breakup data and fix the ANC Co by minimizing 
X 2 relative to the ANC; (iv) determine the error ec in 
Co through xLin + 1> ( v ) calculate the neutron capture 
cross section corresponding to Co ± ec ■ 

There are a few details that can be a source of un- 
certainty to this procedure. As mentioned in the intro- 
duction, there is a nuclear contribution to the breakup 
reaction which is included in our CDCC calculations. 
Fragment-target optical potentials are usually fitted to 
elastic scattering to reduce the ambiguity. In this partic- 
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FIG. 1: Differential cross section with respect to energy for 
15 C -> 14 C+n breakup on 208 Pb. The data are from Ref. U 
and the lines represent the cross sections obtained from each 
of the potential sets in Table [T] 
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FIG. 2: Capture cross sections, multiplied by the energy fac- 
tor E^ 1 / 2 , versus neutron energy. The shaded area are cross 
sections obtained from the RIKEN data jjj] and the black cir- 
cles are the latest direct measurements [3]. 



ular example we use the neutron optical potential from 
Perey and Perey [13] and the 14 C optic al p otential was 
taken from 10 Be+ 208 Pb scattering, as in [23|], with a mass 
scaled radius. It has been shown [24[ that for the breakup 
of 15 C on 208 Pb at 68 MeV/nucleon, the sensitivity to de- 
tails of these potentials is small. 

Finally there is the issue of the single particle structure 
of the projectile. For other cases, the single particle ap- 
proximation has been shown to carry a significant error 
|2o| . As discussed in previous works (see for example fll| 
and Ref. within), the ground state of 15 C is considered 
to be a very good single particle case with spectroscopic 
factor very close to S Sl/2 = 1. Even if this were not the 
case, when both the low energy capture and breakup are 
completely peripheral, the s-wave component dominates 
the cross section. This part is directly proportional to the 
bound state ANC, in such a way that the dependence on 
the single particle parameters and the spectroscopic fac- 
tor is negligible once the ANC is fixed C 2 s — S Sl/2 b?, Sl/2 - 

As mentioned above, in 15 C there is also a d-wave excited 
state. The contribution of this state to the neutron cap- 
ture is small and uncertainty in the structure of this state 
does not affect the errors bars. 

Results. We present in Fig. Q] the results for 
the breakup cross sections of 15 C on 208 Pb at 68 
MeV/nucleon, calculated within CDCC using the model 
space described in the previous section. The shape of the 
distribution compares well with the data. Most impor- 
tantly, the peak of the cross section scales linearly with 
the the ANC-squared. For each of the theoretical curves, 
X 2 was calculated and a quadratic relation with the ANC 
was determined. By minimizing the x 2 i an ANC was 
fixed Co = 1.28 ± 0.01 fm -1 / 2 , the error bar correspond- 
ing to Xmin + 1- These allowed ANC values produce a 
range of possible neutron capture cross sections, shown 
by the shaded area in Fig. [5] 

Plotted in Fig. [5] is the range for <r nn E~ 1 / 2 based on 



the RIKEN Coulomb dissociation data, compared with 
the data from the latest direct measurements 0- The 
agreement is very good. Note that the lowest energy 
point in Fig. [2] at 23 keV does not correspond to a mono- 
energetic neutron measurement. The neutrons at this 
energy have a Maxwellian distribution, so an averaged 
cross section is obtained. For 23 keV too, the prediction 
based on the RIKEN Coulomb dissociation data (7.0±0.2 
//b) compares well with the direct measurement (7.1 ± 
0.5 fib). For the purpose of the comparison in Fig. [2 
we multiplied the 23 keV data point by 0.67, which is 
the factor one obtains assuming a perfect E -1 / 2 energy 
dependence in the cross section (valid at this low energy) . 

A lower energy breakup measurement is also available 
Q . A direct comparison of theoretical cross sections with 
the experimental data was not possible for this experi- 
ment due to a non-linear energy response function of the 
detectors. Therefore the theoretical cross sections had to 
be folded with the detector efficiency in order to compare 
with the data. The analysis in Ref. [8| suggested an (n, 7) 
cross section approximately half that found in the analy- 
sis on the RIKEN data presented in the previous section 
and other direct and indirect measurements p, [lj| UU • 
Here we present CDCC calculations which test the as- 
sumptions that appeared in the analysis of Ref. [§| . For 
the pur pos e of this study we use the single particle model 
of Ref. [lj| , with again the Perey-Perey [22| neutron-Pb 
potential and the same optical potential for the core as 
in the previous section. 

The first important assumption in Ref. [8[ is that the 
nuclear contribution can be subtracted from the data to 
leave a Coulomb only cross section. This was attempted 
by measuring the breakup data on a range of targets 
from the heavy Pb down to the light C target. Assump- 
tions were made on how the Coulomb and nuclear cross 
sections scale with target mass, and by adding them in- 
coherently, a least squares fit of the data was performed 
to estimate the relative cross sections so that the nuclear 
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FIG. 3: (Color online) Nuclear and Coulomb interference in 
the Coulomb breakup of 15 C at 35 MeV/nucleon. 
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FIG. 4: (Color online) Partial wave decomposition of the cross 
section in the Coulomb breakup of 15 C at 35 MeV/nucleon. 

part could be subtracted. 

In Fig. [3] we show the breakup energy distribution for 
the full calculation (solid line), including both nuclear 
and Coulomb, Coulomb only (dotted line), nuclear only 
(dashed line) and the incoherent sum of Coulomb and 
nuclear (dot-dashed line). The nuclear contribution is 
not negligible and interference effects are large, in agree- 
ment with the result of [I?) • Most importantly, the shape 
of the distribution is changed when interference is taken 
into account. At low energies the full calculation includ- 
ing nuclear and Coulomb coherently is actually less than 
the Coulomb only calculation. 

The other main assumptions appear in the analysis of 
the detector efficiencies. In order to calculate the effi- 
ciencies, a cross section must be entered into the Monte 
Carlo simulation of the detector response. The cross sec- 
tion used in Ref. @ was obtained from first order pertur- 
bation theory. At this lower energy, the main assump- 
tions of this theory over simplify the reaction mechanism, 
namely that the breakup proceeds via a single step El 
transition. 

In our CDCC calculations we have included all multi- 
polarities up to Q = 3 and 14 C+n partial waves up to 
I < 3. In Fig.0]we show the contribution of all the partial 
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FIG. 5: (Color online) Higher-order effects in the Coulomb 
breakup of 15 C at 35 MeV/nucleon. 

waves. Although the p-wave is dominant (and is mainly 
El but would contain some nuclear contributions too), 
one should not neglect the other partial waves as it af- 
fects the normalization and the shape of the distribution. 

In Fig. [5] we compare the full CDCC calculations with 
the corresponding one-step calculation. This corresponds 
to a DWBA calculation where the optical potentials used 
are the CDCC coupling potentials as discussed in [n| ■ It 
had already been shown [l6| that multistep effects can be 
very large at energies close to the Coulomb barrier. Re- 
sults presented in Fig. [5] show that at 35 MeV/nucleon 
multistep processes are also important and affect the nor- 
malization of the cross section. 

These assumptions together overestimate the cross sec- 
tion used in the detector simulations, which in turn un- 
derestimate the efficiencies. Hence the theoretical cross 
sections which fit the data, once folded with the detec- 
tor response functions come out much smaller than they 
should. This goes a long way to explaining why the cap- 
ture cross sections derived from this data are approxi- 
mately half that seen in other experiments. This outlines 
the need for good theoretical-experimental communica- 
tion beyond the usual comparison at cross section level. 

To summarize, in the analysis of Q it was assumed 
that: (i) nuclear-Coulomb interference was insignificant; 
(ii) the outgoing neutrons were all in p-waves and (iii) 
the breakup process occurred in one-step. Given that 
these three assumptions are not correct, a re-analysis of 
this data would be highly desirable. 

Conclusions. A procedure for determining the direct 
neutron capture reaction from intermediate Coulomb dis- 
sociation data is presented. The procedure is valid when- 
ever the neutron capture is fully peripheral. Full CDCC 
calculations for the breakup process are compared to 
Coulomb dissociation data and a range of allowed asymp- 
totic normalization coefficients is extracted from x 2 mini- 
mization. Neutron capture cross sections consistent with 
this range of ANCs are determined. 

The method is applied to 14 C(n, 7 ) 15 C. The Coulomb 
dissociation data of Q is analyzed with CDCC and the 
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ANC C = 1.28 ± 0.01 fin" 1 / 2 is obtained for the 15 C 
ground state. We show that the corresponding (n, 7) 
cross sections are consistent with direct measurements 
, providing comparable accurac y, a nd consistent with 
the results from mirror symmetry [ll| . Previous discrep- 
ancies between all three methods where shown to be dras- 
tic in Ref. [ll[. These discrepancies have been partially 
resolved by an improved analysis of the direct measure- 
ment by Reifarth et al. 0] , with the exception of Q . Our 
calculations imply that assumptions made in the analysis 
of Horvath et al. @ are not valid. In order to include 
this measurement in the extraction of the 14 C(n,7) 15 C 



cross section, a reanalysis of the data is needed. 
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